Ligand field parameters and the ground state of Fe(II) phthalocyanine 
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A judicious analysis of previously published experimental data leads one to conclude that the 
ground state of iron(II) phthalocyanine is an orbitally degenerate spin triplet a\ g e\^02 ( 3 E g ). The 
ligand field parameters, in relation to Racah's C, are approximately as follows: B20/C = 0.84, 
B40/C = 0.0074. The uniqueness of this result is demonstrated by means of a special diagram in 
the B20/C — B40/C plane (under additional conditions that B44/B40 = 35/3 and B/C — 0.227). 
The system is in a strong-ligand-field regime, which enables the use of single-determinant techniques 
corrected for correlations within the 3d shell of Fe. 

PACS numbers: 



I. INTRODUCTION 

Over several decades the interest in iron (II) phthalo- 
cyanine (FePc) has been motivated by various applica- 
tions as well as by its proximity to iron proteins. In more 
recent times FePc has become popular as a model sys- 
tem for core-level spectroscopy (XAS, XMCD) studies*^— 
Certain progress has been made in creating artificial or- 
dered structures of FePc4^ Perhaps the most significant 
finding of the recent decade was the discovery of an un- 
quenched orbital moment of iron in FePc by means of 
in-field Mossbauer spectroscopy 6 and XMCD. 3 

This has in turn brought about a surge of computa- 
tional activity on FePc. Along with various density- 
functional calculations,— ~— it is worth to mention mul- 
tiplet structure calculations 1,2,16 based on a phenomeno- 
logical model. The main ingredients of the latter ap- 
proach are the Coulomb repulsion, allowed for by way 
of the Slater-Condon parameters, and the crystal (lig- 
and) field (CF) on Fe 2+ . The multiplet calculations^ 2 -^ 
were mainly aimed at simulating x-ray absorption spec- 
tra; however, they produced an interesting by-product. 
This is a map of ground states of Fe 2+ in CF parameter 
space (Fig. 2 of Ref . 1) . An early version of such a dia- 
gram for the point group D^h was produced by Konig and 
Schnakig,— but the idea itself goes back to the classical 
work of Tanabe and Sugano^ who dealt with the cubic 
symmetry. Unfortunately, in the case of D^h one can- 
not plot but 2-dimensional sections of the 3-dimensional 
space of CF parameters, the choice of these sections in 
Refs. [13 and HI being rather suboptimal. Besides, the di- 
agrams in Refs. [I?] and [l| have the disadvantage that CF 
parameters in energy units are plotted on the axes, and 
so the diagrams depend on the Slater-Condon (or Racah) 
parameters employed. As against that, the original work 
of Tanabe and Suganoi^ presented the result in terms 
of a dimensionless ratio of the CF parameter to Racah's 
B, which led to the celebrated series of universal dia- 
grams. Still, Micdcma's diagrams are of interest. They 
have an enigmatic cornered shape, the domain bound- 
aries are piecewise-linear, with repeatedly encountered, 
characteristic slopes. These features of the diagrams have 



so far remained unexplained. 

As regards agreement with experiment, the calcula- 
tions leave much to be desired. Density-functional calcu- 
lations make inconclusive predictions of the ground state. 
Thus, Reynolds and Figgis^ could not decide between 
3 E g and 3 B 2g because the two lie too close in energy. 
Marom and Kronik 9 found either 3 B 2g or 3 Ai g , depend- 
ing on computational details. More recently, Nakamura 
et alj^ 2 - found 3 A\ g in an isolated FePc molecule, but 
3 Eg in a columnar stack of such molecules. Establish- 
ing the symmetry of the ground state does not settle the 
dispute: within the correct 3 E g one should further distin- 
guish between the configurations b\ g e 3 a\ g , as conjectured 
by Dale et al.,— and a\ g e 3 b\ g , found in the multiplet 
calculations, 1,2 The two ground-state configurations lead 
to distinct types of magnetic behavior. 

This work aims at determining the CF parameters 
of Fe(II) phthalocyanine. As we will show below, the 
known experimental facts on that compound (obtained 
by magnetic and spectroscopic measurements) in connec- 
tion with our CF analysis let no choice: there is only one 
domain in the field of CF parameters yielding a ground 
state that does not contradict established knowledge. In 
such a way our calculations resolve the confusing puz- 
zle about the groundstate of FePc that existed for many 
years. As a byproduct we explain the peculiar shape of 
Miedema's diagrams. 

In the following, we consider the 3d 6 configuration in 
a CF of symmetry D^h and allow for Coulomb repulsion 
between the 3d electrons. The CF is a priori assumed 
neither strong nor weak as compared with the Coulomb 
interaction. Hybridization of the Fe 3d orbitals to neigh- 
boring ligands is thought to be included into the relevant 
CF parameters. In that sense it is better to call our the- 
ory ligand field theory instead. But the ligand p orbitals 
are not treated in an explicit way and we keep the term 
CF theory for simplicity. The spin-orbit coupling is ne- 
glected at first (since it is much weaker than either the 
CF or the Coulomb repulsion) but taken into account in 
a later discussion of magnetic properties. 

This paper is organized as follows. In the next sec- 
tion we briefly review the experimental facts that bear 
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on our knowledge of the ground state of FePc and reit- 
erate the current status of this knowledge. Further, in 
Section III, a diagram of ground states of FePc is con- 
structed from numerical calculations. Our diagram is 
similar to that of Ref. [l|, the main two differences being 
that (i) dimensionless coordinates of the Tanabe-Sugano 
type are used and (ii) the section of the 3-dimensional 
space of CF parameters is chosen on the principle that 
the coordination polyhedron is a plane square. In Section 
IV, the same diagram is reproduced analytically, which 
includes explicit expressions for all domain boundaries. 
The subsequent discussion in Section V hinges upon the 
good agreement of the exact (numerical) and approxi- 
mate (analytical) diagrams. The piecewise-linear shape 
of the domain boundaries finds a natural explanation in 
the linearity of the underlying equations. A conclusion is 
made that FePc is in a strong CF mode and approximate 
values of the CF parameters are given (or rather, ratios 
of CF parameters to Racah's C) . The ground-state con- 
figuration turns out to be a\ g e 3 b\ g ( 3 E g ), as in Refs. [Hi. 
Section VI recapitulates the conclusions. 



II. EXPERIMENTAL FACTS AND THEIR 
IMPLICATIONS 

A. Magnetic susceptibility 

As early sources of our knowledge of the ground state 
of FePc one usually cites magnetic susceptibility stud- 
ies of /3-FePc powder— and single crystals^ The ex- 
perimental data of both papers are in reasonable agree- 
ment with each other. At temperatures between 100 and 
300 K the susceptibility follows the Curie- Weiss law with 
fj, cf f ss 3.8 /ie (for powder). This is between the spin-only 
values of fi c s for S = 1 and 5 = 2 (2\/2 (ib ~ 2.8 /ub 
and 2\/6 hb ~ 4.9 /ib, respectively). Below about 20 
K the susceptibility of /3-FePc becomes temperature- 
independent. 

These facts found an explanation in a simple model 
with S — 1 and effective ^-factors employed in both 
works . 20 i 21 The spectrum of the model consists of a sin- 
glet ground state with M$ = and an excited doublet 
with Ms — ±1 situated at ~ 70 cm -1 . It is unclear why 
Dale et ali 2 ^ thought to justify this model by proposing 
H g e l a lg ( 3 Eg) as the ground configuration (and calling 
it an orbital singlet). Their work contains no experi- 
mental evidence of 3 E g being the ground state of FePc. 
Barraclough et ali 2 ^ noticed the discrepancy between the 
orbitally degenerate 3 E g and Dale's assertion that the 
ground state should be an orbital singlet, and postulated 
3 i?2p instead. As pointed out in Ref. [22|, this was no 
proof, 3 Ai g or 3 A 2g could have done equally well. 

The model used in Refs. H3 and [5l| is not without 
its difficulties. So, it cannot explain the presence of an 
excited state (or states) at ~ 300 cm -1 , as pointed out 
in Ref. [5(3. The existence of such an excited state follows 
from the fact that the susceptibility deviates from the 



Curie- Weiss law above room temperature, as observed by 
Lever £2- (A slight downward curvature is also visible in 
X" 1 vs T data obtained more recently on a-FePc£) This 
can be viewed as an argument in favor of 3 E g rather than 
an orbital singlet. The six- fold degenerate 3 E g would be 
split by the spin-orbit interaction into 4 singlets and a 
doublet, the overall splitting being ~ Q ~ 400 cm -1 . The 
observed susceptibility behavior would find a plausible 
explanation if one of the singlets was the ground state, 
the doublet (or a quasi-doublet) was situated at ~ 70 
cm -1 , and a further state (or states) at ~ 300 cm -1 . 

Another difficulty of Dale's triplet model consists in 
the values of the g-factors, which differ significantly from 
2. Thus, Dale et al3 obtain g x = 2.86 (and g\\ = 1.93). 
That is, nearly one Bohr magneton has to come from an 
orbital moment. Such a large orbital contribution is ex- 
plained more naturally by the presence of an unquenched 
orbital moment (i.e., by orbital degeneracy of the ground 
state) rather than by mixing in of excited states. We note 
that Barraclough et al.jS! who assert most emphatically 
the equivalence of their approach to that of Ref. [2(J ob- 
tained an isotropic g-factor, g± = g\\ = 2.64. Generally 
speaking, Barraclough's (/-factors should be more trust- 
worthy, since they were deduced from data measured on 
a single crystal. 21 The difficulty, however, is that, accord- 
ing to Eq. (4) of Ref. [2(| the zero-field splitting must 
vanish for g± = g\\- At the same time, it is emphasized 
that this splitting, ~ 70 cm' 1 , is very large<21 

In any case, it should be regarded as firmly established 
that the susceptibility is a maximum in the plane of the 
FePc molecule^ This conclusion has been recently con- 
firmed in an independent experiment £ As regards the 
ground states conjectured to explain the susceptibility 
data, they cannot be viewed as deduced from experiment. 

B. Other techniques 

An x-ray diffraction experiment of Coppens et alJ^ 
found the occupation numbers of the Fe 3d orbitals in 
FePc: &2g 65e2 .' 13a i g 88 fri ff 75 - On account of covalency, these 
numbers sum up to 5.41 rather than 6. Restoring the nor- 
malization to 6, one has b\ g 3 e 2 ' 36 'a?g 98 6?j 83 . Coppens et 
al. regarded their result as a direct confirmation of Dale's 
conjecture, b\ e 3 a\ { 3 E g ). Yet, the analysis in Ref. 

was limited to spin-triplet states. An unprejudiced 
look at the quintet states, in particular at b% g e^a\ b[ g 
( 5 B2g), suggests a higher degree of agreement with Cop- 
pens' results. However, 5 -629 can be ruled out because 
it would have resulted in too high a magnetic moment, 

jtZeff = 4.9 

Turning now to the optical absorption experiments of 
Stillman and Thomson^ we note that they were carried 
out on FePc solution in dichlorobenzene. This system 
is chemically distinct from either the free FePc molecule 
or a or /3 FePc. Therefore, without casting doubt upon 
Stillman and Thomson's assertion of a 3 ^2g ground state, 
we state merely that their result is not relevant to the 
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system under consideration herein. 

A Mdssbauer spectroscopy study of Filoti et al£ found 
in a-FePc a very large (66 T) hyperfme field on 57 Fe. Un- 
like the usual Fermi's contact field, the hyperfme field in 
a-FePc has a positive sign (meaning H^f ft A*Fe) and can 
only originate from a large unquenched orbital moment. 
The latter was estimated to be about 1^ but no definite 
information about its orientation could be obtained. 

A more recent XMCD experiment of Bartolome et al£ 
found in FePc an orbital moment of 0.53 (i^ lying i n the 
plane of the molecule. In the same work 3 it was demon- 
strated by direct measurements that the plane of the 
molecule contains the easy magnetization direction, in 
agreement with the early finding of Barraclough et alj^i 

To summarize the section, there is no experimental 
evidence of the ground state of FePc being either 3 i?2g or 
3 Ai g or 3 A2 g . Nor do Coppens' data^ provide sufficient 
confirmation for Dale's conjecture of b% g e^a^ g 
All one can say at this point is that it should be a 3 E g 
state endowed with magnetic anisotropy of an easy-plane 
kind. 



Thus, the nearest environment of the iron atom consists 
of four nitrogen atoms making a plane square, the Fe-N 
bonds being aligned with either the x or the y axis. This 
fact enables us to reduce the number of independent CF 
parameters by one. A rather general CF model known as 
the superposition model (see Ref. HH for a comprehen- 
sive review), relates pairs of CF parameters B nm with 
equal n on the basis of shape of the coordination polyhe- 
dron. Omitting the rather straightforward calculations, 
we state the result: for a plane square the superposition 
model demands that 
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B. Hamiltonian matrix 



(4) 



Our calculations dealt with a Hamiltonian consisting 
of the CF ^ and the Coulomb repulsion and operating 
on the 3d 6 configuration. The basis states were taken in 
the form of simple products of one-electron d orbitals, 



III. NUMERICAL CALCULATIONS 

A. Crystal field Hamiltonian 

The CF Hamiltonian operating on a single 3d electron 
in a tetragonal {D^h) environment is written as follows: 



7-Lcf = B20O2 + -B40O4 + S44O4 



(1) 



Here are Stevens' operator equivalents^ in the t- 
representation (I = 2): 0% = 3^-6, 0% = 35^-155^ + 
72, 0\ = + £-); B nm are CF parameters. In older 
literature one sometimes comes across Ballhausen's CF 
parameters.— These are related to the B nm 's in a simple 
way: 



Dq 



Ds = 3B 20 , Dt = %Bu-l2B 40 (2) 



It is well known that the five real d orbitals belong to dis- 
tinct irreducible representations of the point group D^. 
Therefore, in the basis of those orbitals the CF Hamilto- 
nian ([T]) takes a diagonal form, the eigenvalues being^ 



E(d xy ) 

(ydxz,yz ) 

E(d z 2) 
E(d x 2_ y 2) 



E(b2 9 ) 
E(e g ) 
E{a lg ) 
E(h g ) 



6B20 + 12£>4o — I2B44 
— 3i?2o — 48£>4o 
-6B20 + 72£?4o 
6B 20 + 12B 40 + 12B 44 



(3) 



(5) 



i=l 



with rrn = 0, ±1, ±2, and a, = ±1/2. There are ( 6 U ) = 
210 such states in total. 

Nonzero matrix elements of %cf © are of two kinds. 
First of all, there are diagonal matrix elements, given by 



B 



20 



3m,- — 6) 



B 4 



35mf - Ibbmj + 72) (6) 



i=i 



Secondly, there are nonzero matrix elements between the 
states ([S]) that differ in one pair of quantum numbers rrii, 
rrii being —2 in one of the states and +2 in the other one. 
All such matrix elements equal I2S44. 

The matrix elements of the Coulomb repulsion have 
been treated extensively in the literature. Here we follow 
Griffith's fundamental treatise^ Again, there are two 
distinct kinds of nonzero matrix elements. The diagonal 
ones are given by 



E 

fc=0,2,4 



i>3 



(7) 



where F k are the Slater-Condon parameters (k = 0, 2, 4) 
and 



Note that Ballhausen's original equations (5-14) and (5- 
15) need to be augmented with the cubic terms, +6Dq 
and —ADq, respectively, before being converted to the 
Stevens notation by means of Eqs. @. 

So far no restrictions have been imposed on the CF, 
except that it should be compatible with the point group 
D^h- Yet, much more is known about the structure of 
the FePc molecule than just the symmetry of the Fe site. 



4tt 



2fc + 1 



2m 2m' 1 k,m—m 



(8) 



The integral in Eq. is known as the Gaunt coefficient. 
Numerical values of c^ m / were taken from Table 4.4 of 
Griffith's book,— The inner sum in Eq. ([7]) is taken over 
all 15 pairs of filled d orbitals. The first term in brackets 
describes the so-called Coulomb contribution, while the 
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second one, relevant to pairs of orbitals with parallel spins 
only, is the exchange contribution. 

The Coulomb repulsion also has off-diagonal matrix 
elements. These are nonzero only between the states 
with equal Ml and Ms that differ in two occupied d or- 
bitals, say, \mnan} and \mjiUj\) in State #1, as against 
\mi2<Ji2) and \rrij2crj2) in State #2. It must hold that 
mn+rriji = mi2+mj2 and <Jn+(Tji = <Ji2 + ®j2- The ma- 
trix element between the states #1 and #2 is expressed 
as follows 



E < 

fc=0,2,4 



c k c k 

1 mantis THjz'mji 



0.1 

0.08 : 
0.06 - 
0.04 - 
0.02 - 
- 
-0.02- 
-0.04- 

-0.08- 
-0.1- 





e,„/c 



Thus, the matrix elements of the Hamiltonian are linear 
combinations of the CF parameters, i?20i B±q, and -B44, 
as well as the Slater-Condon parameters, F°, F 2 , and 
F 4 . The latter are conveniently replaced by the Racah 
parameters, 

F° = A + \C, F 2 = 495 + 7C, F A = f C (10) 

The parameter A is hereafter set to zero, because its only 
effect is to shift the energies of all the states of d n by the 
same amount, An(n — l)/2. 

C. Degeneracy diagram 

The calculation consisted in setting and numerically 
diagonalizing the Hamiltonian matrix for given sets of 
parameters B20, -B40, -B44, B, C, and subsequently de- 
termining the degeneracy of the ground state. Five char- 
acteristic values of degeneracy were encountered: 



1 : 


l A 


S = 


0, no orbital degeneracy 


3 : 


3 A 


S = 


1, no orbital degeneracy 


5 : 


b A 


S = 


2, no orbital degeneracy 


6 : 


3 E 


S = 


1, double orbital degeneracy 


10 : 


5 E 


S = 


2, double orbital degeneracy 



At this stage the ground states are labeled tentatively. 
So A can be anything of the following: A± g , A2 g , -B lg , or 
B2 g , which we are unable to distinguish. On the other 
hand, 3 E= 3 E g and 5 E= 5 E g as will be explained in Sec- 
tion IV. 

The construction of the diagram (Figure [lj was or- 
ganized as follows. All energies were expressed in the 
units of the Racah parameter C. The ratios B20/C and 
B40/C were treated as independent variables defined on 
a dense mesh. In the spirit of Tanabe and Sugano, 18 
the ratio B/C was fixed to a value appropriate for Fe 2+ , 
B/C = 0.227, as in Table 7.3 of Ref. |1 The CF pa- 
rameter B44 was not regarded as an independent one. 
Rather, it was found from Eq. (U), as prescribed by the 
superposition models As a result, the B20/C — B^/C 



FIG. 1: Partition of CF parameter space among differently 
degenerate ground states, as found from numerical calcula- 
tions in the absence of spin-orbit interaction. The possible 
ground states are denoted according to their total spin and 
the absence ( 2S+1 A) or presence ( 2S+1 E) of orbital degener- 
acy. They will be further specified in Figure [3] 

plane was partitioned into domains of five different kinds, 
according to the degeneracy found at each point. The 
diagram (Figure [1]) has a cornered shape reminiscent of 
the diagrams in Refs. [Hand [13. The domain boundaries 
appear straight lines with characteristic slopes. Several 
sets of parallel lines are encountered. The central part of 
the diagram is an area of weak CF; in compliance with 
Hund's first rule, the ground state here has S = 2. The 
periphery of Figure[T]is a region of strong CF; here 5 = 
or 1. 

The numerical calculations have the advantage of pro- 
ducing an immediate graphical result. However, it is not 
easy to analyze the character of a ground state expressed 
in a 210-dimensional basis. Of special interest to us is 
3 'E g , which appears in six non-adjacent domains in Fig- 
ure [TJ So we would like to know if those 3 E g are similar 
or distinct. Furthermore, we would like to find out the 
origin of the cornered shape of the diagram, why the 
boundaries are straight and the slopes repeated. Finally, 
we pose a question, how the diagram would change if 
B/C and/or B44/B40 were different to the ones used so 
far. Answers to the above questions should be sought by 
means of analytical calculations. 



IV. ANALYTICAL TREATMENT 

A. Weak crystal field: quintet states 

In the weak-field approximation the CF is treated as a 
perturbation with respect to the intra-atomic Coulomb 
interaction, whose eigenstates are spectral terms with 
certain L and S. Since the CF acts on spatial but not on 
spin variables, terms with different S do not mix together 
(as long as the spin-orbit coupling is neglected). In the 
d 6 configuration there is a single quintet term, 5 D, whose 
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Coulomb energy is 27 

^Coulomb = -35S + 7C (11) 

The remaining task consists in diagonalizing the CF 
Hamiltonian (JXJ) on the wave functions of 5 D, since there 
are no other terms with S = 2. To this end, it is conve- 
nient to interpret Eq. (JT|) in a slightly different way than 
it was done in Section III. Namely, O™ are now regarded 
as Stevens' operators in the L representation (L = 2): 
2 = 2>L 2 Z — 6 etc. Since 5 D contains a single d electron 
above a closed semi-shell, it is only this one electron that 
is exposed to the CF. Therefore, L — I and the coeffi- 
cients B nm in Eq. (JT]) are the same in both representa- 
tions. So we can simply take over the one-electron CF 
energies ©. In doing so, we capitalize the irrep labels, 
to indicate that they now refer to many-electron states, 
and indicate the multiplicity 5. We also prefix -©coulomb 
(ITTj) . The resulting energies of the quintet states are as 
follows: 

E( 5 B 2g ) = -355 + 7C + 6B 20 + l2B ia - 12B 44 (Ql) 
E( 5 E g ) = -35B + 7C- W 2Q - 48S 40 (Q2) 
E( 5 A lg ) = -355 + 7C - 6B 2Q + 72B 40 (Q3) 
E( 5 B lg ) = -355 + 7C + 6B 2Q + 12B A0 + l2B 4i (Q4) 



fitlj/8,4 




FIG. 2: Partition of the B20/B44 — B40/B44 plane among all 
possible permutations of the four one-electron CF levels. 



B. Strong crystal field: singlet states 

In the strong-CF approximation the zeroth-order 
states are constructed from one-electron eigenstates of 
the CF Hamiltonian (JTJ , then their energies are corrected 
for the Coulomb repulsion. A first question that arises 
is: which six one-electron d states are filled in a CF of 
symmetry D 4 hl To give a possibly general answer, it is 
convenient to express all relevant energies in the units 
of -B44. Thus, the one-electron CF energies fl3J) are di- 
vided by -B44. Then, equating pairs of the so modified 
expressions, one obtains 5 equations linear in B 2 q/B 44 
and -B40/-B44. The corresponding lines in the parameter 
plane B 2 q /B44 — B40 j B44 (Figure [5]) are loci of points 
where the sequence of CF levels changes. For example, 
the levels a\ g and b\ g cross over on a line decribed by 

- B 20 /B 4 4 + 5B40/B44 = 1 (12) 

as readily obtained by equating the last two Eqs. ([3]). 
Equation (fT2j) describes the upper one of the two par- 
allel lines in Figure [2J the lower line arises from the 
condition E(a\ g ) = E(b 2g ). Likewise, the equation 
E( e g) = E(b\ g< 2g) generates a pair of parallel lines with 
a negative slope. A single line passing through the ori- 
gin is produced by the relation E(ai g ) = E(e g ). Finally, 
the equation E(b\ g ) = E(b 2g ) leads to no line; this is 
why there are 5 solid lines in Figure rather than 6 as 
expected combinatorially. The CF levels b\ g and b 2g do 
not swap at any inner point of Figure [51 but do so at 
infinity, where -B44 changes sign. Therefore, b 2 (a short 
for b 2g ) stands always to the left of b\ (bi g ) in the level 



sequences indicated within each one of the 12 domains. 
The sequence labels, read from left to right, name the 
CF levels in order of ascending energy if B44 > 0, and in 
order of descending energy if -B44 < 0. 

Up to this point no restrictions have been placed on 
the CF, apart from those imposed by the Z?4^ symmetry. 
Now we do restrict the CF by demanding that it must 
additionally comply with the superposition model, Eq. 
(j4|). This implies that the system is now bound to the 
horizontal dashed line in Figure [2] The dashed line cuts 
through six domains. The corresponding intervals on the 
abscissa axis are numbered 1 to 6 in order of ascending 
B 2 o/B44, with B44 > 0. Thus, the interval #1 stands 
for B 20 /B 4 4 < #2 means that -|f < B 2 o/-B 44 < 

— j etc. The same intervals, but with B44 < 0, will be 
referred to by overscore numbers 1 to 6. 

Now, examining the CF level sequences in the above 
12 intervals, one encounters only 3 situations where there 
is a CF gap between the highest occupied and the low- 
est unoccupied orbitals. The corresponding ground-state 
electronic configurations are as follows: 

b%e g intervals #2, #3, #4 

a 2 lg e 4 intervals #5, #6, #1 (13) 

a lg b ig b lg intervals #4, #5 

Within the above intervals of B 2 q /B44 the ground state 
is a singlet, provided that the CF is sufficiently strong. 
In all other cases the Fermi level is caught at the par- 
tially occupied quadruply degenerate e g level and there 
is a possibility of both a singlet and a (spin) triplet 
ground state with the same energy. A subsequent al- 
lowance for intra-atomic (Hund's) exchange makes the 
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triplet states energetically more favorable than the sin- 
glet ones. Therefore, singlets are no viable candidates 
for ground state in all intervals where triplets with the 
same CF energy are possible. In such intervals only the 
triplets will be considered (in the next subsection). 

Conversely, in the three cases where there are viable 
singlet states (fT3|) , competing triplet states will be taken 
into consideration as well, constructed from excited CF 
configurations. Such triplets still have a chance of be- 
coming ground state on account of Hund's exchange in 
situations where the CF is not strong enough, near inter- 
val boundaries, etc. 

Let us now turn to our direct task — computing the 
energies of the singlet states ([T^f. The CF energies are 
computed most readily, by summing up the energies of 
the six occupied one-electron states as given by Eqs. ([3]). 
First-order correlation corrections are then computed fol- 
lowing Slater's prescriptions^ for each pair of occupied d 
states a so-called Coulomb integral J(di,d2) is added; a 
further exchange contribution K{d\,d,2) is deducted for 
pairs with equal spins. J's and K 's between the real d or- 
bitals were expressed in terms of the Racah parameters 
by Griffith, see Table A26 of his booki^l The resulting 
singlet energies are as follows: 

E ( b lg e t) = -I68B40 - 24S 44 - 30B + 15C (SI) 
E{af g e l ) = -24B 20 - 48B 40 + 10B + 15C (S2) 
E (al g b\ g bl g ) = 12B 20 + 192B 40 - 20B + 15C (S3) 

We note that for low-lying single-product states, such as 
those considered in this work, the factor of C depends 
solely on S and is given by 

(factor of C) = (S max + l) 2 - (S + l) 2 (14) 

where S ma x = n/2 is the hypothetical maximum spin of 
n electrons in the absence of the Pauli principle. For 
fewer than six d electrons, states with S = S ma x are 
allowed and their energies have no contribution in C. 
For d 6 , S'max = 3 and the factors of C equal 15, 12, 
and 7 for S = 0, 1, and 2, respectively. Equation (jT4"|) 
is a consequence of the great simplicity acquired by the 
Coulomb and exchange integrals when A and B are set 
to zero: 

J(ck, dj) = K(di, dj) = (1 + 2<%)C, 

cf. Table A26 of Ref. H3- No simple relations are known 
for the factors of B, which have to be calculated in each 
case separately. 

C. Strong crystal field: triplet states 

Construction and finding the energies of the (spin) 
triplet states are carried out in a similar fashion. One pe- 
culiarity is the large number of triplets (9 in total) , which 
have to be constructed for all 12 intervals of B20/B44. 
Where no triplet state is permitted by the ground CF 



configuration, the first excited configuration will be con- 
sidered instead. 

We proceed from the interval #1, B20/B44 < — 12, 
B44 > 0. Here (as well as in the interval #6, B20/B44 > 
B44 < 0) the ground CF configuration is b\ g b\ g e 2 g , 
which allows one triplet state, c? 2 2 _ id 2 xy S xz S yzl as well 
as three singlet ones. According to the first Hund's rule, 
it will be the triplet that will become ground state upon 
allowance for the Coulomb interaction. (It was for this 
reason that the singlets were left out in the previous sub- 
section.) It is easy to figure out that the triplet state 
under consideration has the symmetry of the product xy, 
that is 3 B2 g - Its energy is computed following the same 
prescription as in the previous subsection and equals 

E{b\ g bl g e^) = 18B 20 - 48B 40 - 9B + 12C (Tl) 

Let us move to the interval #2, — 1& < B20/B44 < — |, 
-B44 > 0. The ground CF configuration, b 2 g e g , consists 
of fully occupied orbitals and is necessarily a singlet. To 
construct a spin triplet state, one spin-down electron is 
promoted, with a simultaneous reversal of spin, from the 
e g orbital to the first unoccupied CF level bi g . The result 
is either dl y d 2 xz d\ z d\ 2 ^ or d 2 xy d 2 yz dl z dl 2 _ y2 . This is a 
doubly orbitally degenerate state 3 E g . Its energy is 

E{bl g efH\ g ) = 9£ 2 o-1085 4 o-12B44-24B+12C (T2) 

Proceeding as before, we find that the most favorable 
spin triplet state in the interval #3 is another 3 E g , whose 
energy is given by 

E(b 2 2 g el^al g ) = -35 2 o-48B4o-24B 4 4-28B+12C (T3) 

The remaining six spin-triplet states include: a 3 -B 2 g 
in the intervals #4 and #5, with 

E(epl g a[ g ) = -12B20-108B40-12544-22B+12C ( T4 ) 
a 3 E g in the interval #6, with 

E{a\ g ef^bl g ) = -I5B20+I2B40-I2B44-UB+I2C (T5) 
a 3 E g in the interval #1, with 

E(a 2 lg el l H\ g ) = -155 20 +12B 4 o+125 4 4-14B+12C (T6) 
a 3 B2 g in the intervals #2 and #3, with 
E{a\ g b\ g e^) = -6B 2a +72B 40 +24B44 - 29B+12C (T7) 
a 3 Eg in the interval #4, with 

E{a\ g b\ g bl g el) = 3B 20 +132B 4 o+12B44-295+12C ( T8 ) 

and a 3 E g in the interval #5, with 

E{b\gbl g a\ g el) = 15B 20 + 72B 40 - 13B + 12C (T9) 

Note that the factor of C in Eqs. (T1-T9) is invariably 
12, as follows from Eq. (fT4"|) with S = 1 and S'max = 3. 
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FIG. 3: Partition of CF parameter space among possible 
ground states of FePc. The labels are mnemonic: so Si is 
a singlet whose energy is given by Eq. (SI). The borderlines 
are as described by Eq. (|17p with the coefficients of Table I. 



D. 



The B 2 o/C - B40/C diagram 



The search for the ground state consists in a systematic 
comparison of energies of pairs of candidate states, as 
given by Eqs. (Q1-Q4, S1-S2, T1-T9). For example, 
equating (Ql) to (Tl) results in 



12B 



20 



60S 



40 



12B. 



44 



265 + 5C 



(15) 



Eliminating B44 by means of Eq. Q and dividing the 
result by C, one arrives at an equation of a straight line 
in the plane of the parameters B20/C and -B40/C: 



Bao Bon 

_E = -0.15-^ 



0.0625 - 0.325 



B 

C 



(16) 



Left of this line there should be a domain where the 
ground state is the triplet T\ ( 3 i?2 S or b^gh^e^), right 
of the line, towards the origin, lies the domain where the 
ground state is the quintet Qi ( 5 i?2g)- In the spirit of 
Tanabe-Sugano, the ratio B/C is fixed, B/C — 0.227, as 
in Table 7.3 of Ref. [M 

Proceeding as above, one obtains equations for all 33 
borderlines appearing in the diagram, Figure [3l (It suf- 
fices to consider pairs of states belonging to the same, or 
perhaps, to adjacent intervals of B20/B44.) By analogy 
with Eq. (|16l) . these expressions are presented as 



TABLE I: Values of coefficients in Eq. JT7J). 



label 






y 


T n, 


Of ZU 


1 /I fi 

— i / 1U 


1 ^ 1 AO 

lO/ 'lU 


-t 2^1 


1 /zLO 


1 /9/1 

± / Z4: 


11 /1 9n 

11/ _LZU 


T T. 


q /9nn 

— \) 1 zuu 


n 
u 


^ //in 

— Of 4U 


01 12 


q /9nn 

— of ZUU 


^ /900 
/ zuu 


^ /1 nn 

— Of 1UU 




q /1 fin 
— of 10U 


1 //in 

1 / 4U 


1 / 04 


Si To 


1 /zLO 


1 //LO 


1 /fin 

1 f uu 


T 


0/ ou 


^ /9nn 
0/ zuu 


1 /9^ 

— 1/ zo 


T. D 


q /onn 
— y/ zuu 


1 //in 

1 / 4U 


7 /9nn 
/ /zuu 


T, Do 


n 

u 


1 /^fi 
1 / ou 


1 //in 
1 f L ±\j 


T Do 
-t 4 V2 


q /9nn 

y/ zuu 


i //in 

± / 4:U 


1 ^ /9nn 

10/ zuu 


T Do 


^ / I fin 


1 /fi/L 

± / U4: 


1 % /^9n 

10/ ozu 


J 5 V3 


q /9nn 

£7/ zuu 


1 //in 

1 / 4U 


91 /9nn 

Zl / zuu 


J2 V3 


^ /9n 
— 0/ zu 


1/10 


^ /a 
0/0 


0.7; 


q /an 

y / ou 


^ /an 

0/ ou 


^ /1 n 

Of 1U 


>52-i 4 


^ /^n 
0/ ou 


■3 /9nn 
— 0/ zuu 


— 4/ zo 




q /an 
y/ ou 


n 
u 


n 
u 


DoDo 
V2V3 


1 /AO 

1 f 4:U 


n 

u 


n 

U 


ToT 


q /an 
y/ ou 


n 
u 


9 //in 

— O/ 4U 


TT- 
J 4-1 5 


i //in 
1 / 4u 


n 
u 


1 /1 ^ 
— 1/10 


T-Do 
J 6 W3 


Q /an 
y / ou 


— 1 /I fi 

— 1/ 1U 


91 /an 

z ± / ou 


1 7^3 


u 


1 /t^fi 
— 1/ OU 


q /1 /in 

— 0/ 14U 


T~D , 


^ /^n 
0/ ou 


1 //in 

— If 4U 


■3 /1 nn 

0/ 1UU 


S>2 J 6 


q /9nn 
— y/ zuu 


■3 /9nn 
0/ zuu 


/ zo 


T-Tr 


q /9nn 

y/ zuu 


n 
u 


1 /AO 

O / 'lU 




^ /9n 
— 0/ zu 


n 
u 


n 
u 


T^^J . 
-18^4 


1 //in 
1 / 4u 


1 /9/1 
— 1 / 


1 /9n 
— 1/ zu 


T 7 T 8 


9/80 








S3Q4 


-3/20 


-1/5 


-3/8 


SsTg 


9/80 


3/80 


9/80 




9/80 


1/16 


11/40 


S3T9 


1/40 


-1/40 


7/120 


TiQ 4 


3/50 


1/40 


13/100 


TiT 9 


1/40 





1/30 



These are obtained by means of Eq. ((4|) from the inter- 
val boundaries in Figure [2] so B20/B44 = 24/7 leads to 
B20/B40 = 1/40 etc. One exception that is not on the 
list (fTSl) but is encountered twice in the second column 
of Table I is -3/160. 



-B40 
C 



B20 . , , ,,B 
a— + b + b- 



The numerical factors a, b, and b' are listed in Table I. A 
line is referred to by naming the two domains it separates. 
Remarkably, one finds in the second column of Table I 
repeatedly six characteristic slopes, 



0, 



_9 JL A JL and -L 

200 ' 20 ' 80 ' 50 ' 40 ' 



(17) V. DISCUSSION 

In the preceeding section we constructed a diagram of 
ground states of FePc in the absence of spin-orbit cou- 
pling, Figure [3] A total of 16 distinct ground states are 
present in the diagram: 3 singlets (Si — S3), 9 spin triplets 
(Ti — T g ), and 4 spin quintets (Qi — Q4). The respec- 

(18) tive energies are given by Eqs. (S1-S3, T1-T9, Q1-Q4). 
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Explicit expressions were derived for the domain bound- 
aries, Eq. (|T7|) and Table I. The boundaries are segments 
of straight lines, which is a consequence of the linearity 
of Eqs. (S1-S3, T1-T9, Q1-Q4). This gives the diagram 
its peculiar cornered shape, with characteristic, repeated 
slopes. As clear from the structure of Eq. (fl7|) . the slopes 
do not depend on the ratio B/C. Taking a slightly differ- 
ent B/C would shift the domain boundaries somewhat, 
but will not affect their slopes. As against that, the slopes 
will change if the ratio B^/B^ deviates from the value 
prescribed by the superposition model, Eq. More- 
over, such a deviation of -B44/B40 from 35/3 may lead 
to a loss of parallelity of certain boundary lines. For 
example, from a simple analysis of the one-electron CF 
energies <j3j> one finds 

(slope of TiQi) = (slope of Q 3 Q 4 ) = (5 - B^/B^y 1 
(slope of S 3 Qi) = {2B M /B 4a - 30)" 1 
(slope of S2Q3) = -3/20 

Apparently the above lines are only parallel if the condi- 
tion (Q} is fulfilled. In reality the superposition model is 
an approximation and small deviations from Eq. Q are 
to be expected. In the above example, the borderlines 
T\Q\ and QsQa will remain parallel exactly, while the 
others only approximately. A more extensive analysis of 
this matter is beyond the scope of the present work. 

On the whole, the diagram constructed analytically 
(Figure [3]) is remarkably similar to that calculated nu- 
merically (Figure [IJ . We take it as a sign of validity of 
the strong-CF approximation used to compute the ener- 
gies of the singlet and triplet states, Eqs. (S1-S3, T1-T9). 
(N.B. The quintet energies (Q1-Q4) are essentially exact, 
without relying on the weakness of the CF.) This demon- 
strates the applicability of techniques based on single- 
determinant wave functions, even though it is important 
to allow for correlations (nonzero B and C). 

Our next task is to locate the standpoint of FePc in the 
diagrams of FiguresQ]and[3] In the subsequent discussion 
the domain boundaries are assumed to be positioned as 
in the more accurate Figure [TJ whereas the ground states 
associated with the domains are as constructed analyti- 
cally and indicated in Figure[3J The search can be limited 
to an acute angle adjacent to the abscissa axis, within the 
first quadrant of Figure [3] 

< B40 < 0.45B 20 (19) 

Indeed, the d x 2_ y 2 orbital of Fe overlaps most strongly 
with the ligand orbitals and therefore has a much higher 
energy than the other 3d orbitals, in particular, d xy . By 
Eqs. (0, E{d x i_ y i) - E(d xy ) = 245 44 > 0, whence 
by Eq. ([J}, B40 > 0. To prove the right-hand part of 
the double inequality ([T9|) . one should rewrite the CF 
Hamiltonian (fT]), taken in conjunction with Eq. as a 
classical anisotropy energy, 

E a = B 2Q (3 cos 2 9-1) 
+B 4 o(35 cos 4 0- 30 cos 2 9 + 3 + f sin 4 9 cos 40) 
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FIG. 4: Temperature dependence of reciprocal susceptibility. 
Closed circles are experimental data,— solid line is 1.22 times 
the calculated x _1 - 



and demand that 9 = tt/2, <f> — tt/4 be a local minimum. 
This is to account for the well established fact that the 
easy magnetization direction lies in the plane of the FePc 
molecule^i2i 

A further experimental fact to take into consideration 
is that the ground state is a spin triplet (S — 1) and 
that is is orbitally degenerate ( 3 i? ff ).— £ Within the sector 
defined by the condition (TT9|) there are only two domains 
where 3 E g is the ground state — a quadrangle T3 and a 
triangle T5 . We carried out an extensive numerical study 
of the magnetic susceptibility (with due allowance for the 
spin-orbit coupling) and found that X||(^) > XX (T) ev- 
erywhere within T 3 , but X\\( T ) < Xx(T) inside T 5 . (Here 
the subscript "||" refers to the direction parallel to the 
4-fold symmetry axis.) One has to conclude, therefore, 
that the standpoint of FePc in Figure [3] lies inside the 
triangle T5. The corresponding ground-state configura- 
tion is alge^blg, cf. Eq. (T5). It is distinct from the 
configuration T 3 , or 6| 9 e g^ a i fl ) postulated by Dale et 
al,2£ and adopted by Filoti et al£ On a simple model the 
latter authors have demonstrated that T3 has necessar- 
ily an easy-axis anisotropy, which agrees with our analy- 
sis. The experiment f^^i however, insists on an easy-plane 
anisotropy and so T3 has to be definitively abandoned. 
After all, Dale's choice of T 3 was a mere conjecture, with- 
out a sufficient experimental foundation. It should also 
be noted that Miedema et all 1 -* 2 - proceeded from the cor- 
rect ground-state configuration T5, even though they did 
not explain their choice. 

The difference between the two 3 E g configurations is 
easy to understand. In both cases there is one e g hole; the 
two real e g orbitals (d xz and d yz ) can be combined to give 
states with £ z = ±1. An extra singly occupied orbital in 
T3 is aig (rf z a), with £ z — 0. Therefore, the z-component 
of the total orbital moment is L z = ±1 and the spin- 
orbit coupling leads to an easy-axis anisotropy in T3. In 
T5 it is the d xy orbital (629) that is singly occupied and 
the situation is quite different. Now three orbitals, d xy , 
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d xz , and d yz , are accessible to the holes. If the e g and 
&2 g levels were perfectly degenerate, there would be no 
anisotropy at all. The fact that the degeneracy is lifted 
results in a weak easy-plane anisotropy, such as the one 
observed. 

We undertook an attempt to refine the position of 
the system inside the triangle T5 on the basis of the 
available susceptibility data<21 We find that the most 
likely standpoint is near the left corner of the triangle, at 
B 20 /C = 0.84, B i0 /C = 0.0074. Powder susceptibility 
was calculated as ^x\\ + §XJ-j with B — 917cm -1 and 
C = 4040 cm" 1 (as in Table 7.3 of Ref. HI). The spin- 
orbit coupling constant £ was set to 400 cm" 1 . The so 
computed susceptibility proved higher than the experi- 
mental one and had to be reduced by a factor of 0.8, to 
make both curves match. (Accordingly, in Fig. 2] the ex- 
perimental reciprocal susceptibility 21 is compared with 
the calculated x~ 1 times 1.22.) The reduction factor 0.8 
can be attributed to covalency, neglected in our model. 

Apart from the rescaling, the calculated x _1 (T) does 
agree with the experiment. In our calculation the sextet 
3 E g is split by the spin-orbit interaction. The ground 
state is a singlet and so is the first excited state, situ- 
ated 20 cm" 1 above the ground state. The second excited 
state, at 52 cm" 1 , is a doublet, followed by two singlets, 
at 165 cm" 1 and 225 cm" 1 . It will be recalled that the 
model spectrum of Refs. Ho] and Hi] consisted of a ground 
singlet and an excited doublet at 64 cm" 1 . The most es- 
sential distinction of our spectrum is the presence of an 



excited singlet at 20 cm 1 . A clue to this point might be 
provided by a measurement of the specific heat. 



VI. CONCLUSION 

Published experimental data suggest that FePc has 
an orbitally degenerate ground state with 5 = 1, the 
easy magnetization direction lying in the plane of the 
molecule. There is a single domain in the CF parame- 
ter space where these conditions are met — the triangle 
T5 in Figure [3] The corresponding ground-state config- 
uration is a\ g e 3 b\ g . The standpoint of FePc is situated 
in the left corner of the triangle, about B20/C = 0.84, 
B 40 /C = 0.0074, whereas B 44 is given by Eq. (g]). This 
point lies in a strong-CF region, where the notion of 
single-determinant states has a certain validity. 
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